% by YLX 6/6 
% 基于 wwb findQ程序
% sby 的 findQ 中 M 变量为何？
function Q = findQ(S)
    global I J P PIE
    m_max = max(sum(S, 2)); % sum of dimention 2 (row)
    Q = zeros(I, m_max);
    patient_num_i = sum(S_test(i, :));
    for i = 1:I
        for m = 1:m_max
            if m > patient_num_i
                Q(i, m) = -1;
                % 矩阵宽度选择最大m，多余位置是否需要标记为-1？
            else
                Q(i, m) = elementQ(S, i, m);
            end
        end
    end
